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O ■ ABSTRACT 

' We report high-resolution radio imaging of the recurrent nova RS Ophiuchi 

(RS Oph) during the first month of the 2006 outburst, using the Very Long 
'i> I Baseline Array (VLBA). Observations made on days 20.8 and 26.8 of the outburst 

^ I show a synchrotron-emitting partial shell that is much brighter to the east than to 

the west. Assuming the broad component of the infrared lines corresponds to the 
outermost part of the shell seen by the VLBA, the distance to the source is 2.45 ± 
0.4 kpc. The circular shape and spectral indices of the shell emission challenge 
simple models for the radio structure immediately after the outburst. The second 
epoch also shows an additional, resolved, synchrotron-emitting component well 
to the east of the shell. Its inferred velocity is comparable to the escape speed 



X 

^ ■ from the surface of a high-mass white dwarf. This component was not seen in 

the first epoch. Its appearance may be related to the outflow reaching the edge 
of the nebula created by the red giant wind, which had been re-filling the void 
left by the last outburst in 1985. This eastern component is likely related to the 
jets previously seen in this and other symbiotic stars, and represents the earliest 
clear detection of such a jet, as well as the best case yet for synchrotron emission 
from a white dwarf jet. 

Subject headings: binaries: symbiotic — novae, cataclysmic variables — stars: 
individual (RS Oph) — stars: winds, outflows — radio continuum: stars 
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Introduction 



The most recent major outburst of the symbiotic star and recurre nt nova RS Ophiuch i 
(RS Oph) was discovered near the optical peak on 2006 February 12.83 (JNarumi et aLll2006l ). 
hereafter taken as time zero. RS Oph is a binary system in which a red giant deposit s 
material through its wind onto a white dwarf (WD) companion (see, e.g., lKenyonlll986l ). 
Approximately every 20 years, the system erupts as a nova, as 'runaway' thermonuclear 
burning in the hydrogen-rich white-dwarf envelope leads to a visual brightness increase of 
more than 6 mag in less than 2 days, a sustained period with WD lur ninosity ne ar the 
Eddington limit, and the ejection of material from the WD surface (e.g., iBodd 119871 ). The 
last such eruption was in 1985. The frequency of these outbursts su ggests that the rnass of 



the WD in RS Oph is very close to the Chandrasekhar limit (e.g., iHachisu fc Katol 12001 
Sokoloski et al.ll2006l ). making this a prime candidate for a Type la supernova progenitor 
system. 

Whereas most novae de tected in the r adio emit via thermal free- free bremsstrahlung 
from the shell of ejecta (e.g., ISeaquistlll989l ). post-outburst radio emission from RS Oph in- 
stead emanates primarily from the dense circumstellar material (GSM) that is swept up 
and shocked by the fairly low-mass ejected shell. Observations after the previous out- 
burst in 1985 with the European Very Long Baseline Interferometry Network (EVN), the 
Multi- Element Radio Linked Interferometer Network (MERLIN), and the Very Large Ar- 
ray (VLA) suggested that the radio emission was initially non-thermal at low frequencies, 
with a possible thermal component becoming dominant later in the outburst. The radio 
source co nsisted of a resolved central component plus possib l e bi-polar jets expanding at 1.5 
mas/dav JPadin. Davis. &: Bodel EgSsI: iHiellming et~aDll98d : iTavlor et al.lll989l : lDavislll987l : 
Porcas. Davis, fc GrahamI 119871 ). 



Although jets have been observed from a n increasing nurnber o f accreting white dwarfs 
i n symbiotic stars and s upersoft X-ray sources (jSokoloski et al.ll2004l ) , and even several novae 
( JKawabata et al.ll2006l : llijima fc Essenoglull2003l ). these jets are generally not spatially well 
resolved, and basic questions such as the relationship between jet production and outbursts 
remain. Most jets from white dwarfs have velocities of hundreds to thousands of km s~^. 



and when 



HoUis et al. 



et emission can be spat i ally isolated, it usually appears to be thermal ( e.g. , 
1997l : lKellogg et ahlboOll : iBrocksopp et al.ll2004J : [ciallowav fc Sokoloskill20o3 l . In 



a few cases, however, observations have hinted that particles can be ac celerated sufficiently 



that non-thermal r adio and/or X-ray emission can also be produced (JGrocker et al.l 12001 
Nichols et al]l2007h . 



Here we report the detection of an expanding, synchrotron-emitting ring just weeks after 
the nova explosion, and, between three and four weeks after the outburst, the appearance 
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of a synchrotron-emitting feature to the east of this ring ahgned with the putative jet from 
1985. We describe our observations of the expanding ring and the eastern component in §2. 
We discuss these observations in more detail in 53, and summarize our conclusions in 54. 



2. Observations 

We observed RS Oph with the National Radio Astronomy Observatory'^j (NRAO's) 
Very Long Baseline Array (VLBA) under project code BS167 in two four-hour runs, on 
2006 March 5 and 2006 March 11, corresponding to days 20.8 and 26.8 of the outburst, 
respectively. The first run used all antennas but Hancock, while the second used the entire 
10-antenna array. We observed the source in two frequency bands, consisting of 64 MHz in 
each of the two dual polarizations, centered on 1.6675 and 4.9875 GHz, recorded at a rate 
of 256 megabits per second. To remove atmospheric and instrumental effects, we employed 
phase referencing using the nearby calibrators J1745— 0753 (1?67 away; primary calibrator at 
5.0 GHz) and J1743— 0305 (3?27 away; primary calibrator at 1.7 GHz). The positio n of the 
former was taken from the Second VLBA Calibrator Survey ( Fomalont et al.ll2003l ). with a 
quoted (rms) accuracy of 0.50milliarcseconds; t he position of th e latter was taken from the 
International Celestial Reference Frame (ICRF: iFey et al.ll2004l ). with a quoted accuracy of 
0.77milliarcseconds. We correlated all four polarization pairs using the VLBA correlator in 
Socorro, NM, and calibrated, imaged, and self-calibrated the data in the sta ndard fashion 
using NRAO's Astronomical Image Processing System (AIPS; iGreiseru 120031 ). Phase refer- 
encing aligned the 5.0 GHz images to better than 0.5 mas. The 1.7 GHz images differed by 
~ 8 mas, rather larger than expected from residual ionospheric effects, but not alarmingly 
so, given the poor weather during the first observations. We shifted the lower frequency 
images to align the centers of the rings with those seen at 5.0 GHz. 

Figure [1] shows the total intensity images. To a limit of a few percent (or the rms noise, 
whichever was higher), no linear or circular polarization was detected in either epoch, at 
either frequency. There may be some missing radio emission on scales much larger than 
sampled by the VLBA - the corrected 5 GHz flux density on day 20.8 (~ 32 mJy; see below) 
is somewhat lower t han that seen by MERLIN at 6 GHz on the same day (about 40 mJy; 
O'Brien et aDl2006ah . 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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2.1. The Expanding Ring 



All four images are dominated by a ring which is much brighter to the east than to the 
west. Despite this brightness asymmetry, the ring is consistent with being circular (Fig.[2]). 
In particular, a two-dimensional ring with an inclinat ion corresponding to that of the binary 



orbit (~ 50°: iBrandi et al.l 120071 : iQuiroga et al.ll2003l ) provides a significantly worse fit than 



a face-on ring or sphere. These results are based on simulated VLBA observations of circu- 
lar and elliptical rings of various thicknesses and axis ratios, with brightness distributions 
around the ring similar to that observed. These simulations take account of such important 
interferometric uncertainties as the effects of limited uv-coverage and non-linear processing, 
including deconvolution in particular. 

The ring center (as measured at 5 GHz on both days) at ( J2000) is: 

17''50'"13!1583, -06°42'28'.'517, 

with an uncertainty of at most one milliarcsecond (mas). We checked for systematic errors in 
our astrometry by phase referencing our secondary calibrator (J1743— 0350) to our primary 
calibrator (J1745— 0753) at 5 GHz. The resulting position agreed with that given in the 
ICRF to 0.4 mas in right ascension and 0.6 mas in declination - well within the la error bars 
reported in the ICRF. Since these calibrators are ~ 4° apart, while RS Oph is separated 
from the primary calibrator by only 1?67, systematic errors are far below other uncertainties 
in the measurement. Our VLBI position for RS Oph is however ~ 100 mas from the optical 



posit ion from the Naval Observatory Merged Astrometric Dataset (NOMAD: IZacharias et al. 
2004h : 17''50"'13!1628 ± 9mas, -06°42'28'.'496 ± 11 mas. This optical position corrects for 



proper motion to the epoch 2006.19, and is shown as a cross in Fig. [H There is no obvious 
explanation for the radio/optical offset. 

The width of the ring is barely resolved in the first epoch, and clearly resolved in the 
second, suggesting a growing thickness with time. We fit the observed ring as a circle with 
an asymmetric brightness distribution, thickened by convolution with a symmetric Gaussian, 
and then convolved with the beam (point spread function). The fuU-widths-at-half-maxima 
(FWHMs) for the smoothing gaussians were 2.5^25 '^^s *-"^ ^^^ ^^^^ ^^y^ ^^^ ^-^ ^ 0.5 mas 
on the second. The corresponding radii were 16.25 ± 0.5 mas and 17.5 ± 0.5 mas. These radii 
refer to the mid-points of the annuli; the outermost parts of the rings (mid-point -|- FWHM, 
from the fits) are at ~ 17.5 and ~ 19.7 mas. The error bars quoted here are more nearly 
3cr than la errors, insofar as this has meaning for non-Gaussian error distributions, and the 
difference in widths is thus likely to be real. 

The asymmetry of the ring is much more pronounced at 5 GHz than at 1.7 GHz. This 
difference is not an instrumental effect: simulating 5 GHz observations of the 1.7 GHz model 
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shows that the 5 GHz observations would have shown the western emission clearly, had it 
been as (relatively) strong as at 1.7 GHz. Thus the spectrum of the eastern part of the ring 
is intrinsically relatively flat, while that of the western part falls sharply with frequency. 
Based on the simulations, we multiplied the observed flux density at 5 GHz by a factor 1.07 
on day 20.8, and 1.10 on day 26.8, as a first-order correction for the differing uv-coverage 
between 5 and 1.7GHzo The resulting flux densities at 1.7 and 5 GHz are given in Tabled] 



2.2. The Eastern Component 

Our second set of observations revealed an additional component ~ 52 mas to the east 
of the center of the ring. This new eastern component was present at both 1.7 and 5.0 GHz. 
The 1.7 GHz emission at this location is clearly resolved in the north/south direction, with a 
Gaussian fit yielding a rough deconvolved size of 23.7±2.2mas (58.0 ±5.4 AU at 2.45 kpc [see 
§3.1.2l for a derivation of the distance]). As with the ring, we simulated 5.0 GHz observations 
of the 1.7 GHz model, and found that, while the eastern component would appear much 
smaller at 5.0 GHz (as observed), the source spectrum must be intrinsically quite steep (see 
Table [I]). 

This component was not seen in our first set of observations. The 1.7 GHz map shows no 
source to the east of the ring: the 5a flux density upper limit is 0.6mJybeam~^, compared 
to the peak flux density on March 11 of l.SmJybeam^^. The 3cr 5 GHz upper limit is 
approximately 1 mJy beam~^. Further VLBA observations on 2006 Mar 18 through 29 clearly 
detect the eastern component; we postpone discussion of those data to a future paper. 



3. Discussion 
3.1. The Ring 

3.1.1. Expansion Rate and Shape 

The expanding ring is most naturally associated wit h the shock wave resu l ting from the 
outbu r st, and seen i n the s trong early X-ray emission (e.g.. lSokoloski et al.ll2006l : IOsborne et al 
2006al: ISode et al.l l2006b ) and b road optical, infrared, and X-ray lines (e.g., iBuill l2006l : 



Ness et al. 



2006 



Das et al.ll2006bl ). Those data suggest substantial deceleration, as expected 



for a shock sweeping up more and more of the red giant's stellar wind. Figure[3] compares 



^Figures[T]and[2]do not reflect this correction, as it is very likely to vary with position. 
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the radius thus predicted from the optical, IR, and X-ray velocities with that observed in 
the radio. The solid circles represent outer radius estimates from the current data, derived 
through simulatioi is of thick, asy i nmetri c rings as discussed in §2.1[ The open triangles rep- 
resent the radii of lO'Brien et al.l (l2006d ) , measured as half the maximum distance betw een 
peaks in a series of north-south slices. The nominal error bars of lO'Brien et al.l (J2006d ) are 
0.2 mas, but, as they point out, these nominal errors almost certainly underestimate the true 
ones, and we plot Imas error bars in the figure. The systematic difference between O'Brien 
et al.'s measurements and ours, with the former lying some 20% low, is likely due to the dif- 
ferences in fitting procedures, combine d with an evolving surface brightness distribution and 
changing angular resolution (cf. §5.6 of lBartel et al.ll2002l ). Nonetheless, the overall trend of 
radio radius with time, as measured by either group, matches the X-ray/optical/infrared pre- 
dictions nicely. The agreement between the size of the radio ring and the expected location 
of the blastwave also support the conclusion that, unlike the IR continuum e mission, which 



primarily emanates fr om a region within a few mas of the central binary (JMonnier et al. 



20061 : iLane et al.l 120071 ). the radio emission is associated with the expanding shock. 



As shown in Fig.Ej the ring is basically circular. This geometr y is consistent wit h a 
spherical shock projected onto the plane of the sky (see, for instance, iBartel et al.l 120021 . for 
similar observations of supernova 1993J) or a two-dimensional shock in the plane of the sky. 
There is no evidence in the ring shape either for an asymmetric explosion or for any special 
geometry of the circumstellar material. In particular, we s ee no sign that the ring is inclined 



either parallel or perpendicular to the binary orbit (~ 50°: iBrandi et al.ll2007l : iQuiroga et al. 



20031 ). The fractional width of the ring (thickness divided by radius) increases from 0.15" 



-0.06 
-0.15 



on day 20.8, to 0.26 ± 0.03 on day 26.8 (~ 3cr error bars). This breadth suggests that the 
emission arises in the forward rather than the reverse shock, since the latter is expected to 
be much thinner. The change in relative width shows t hat the evolutio n is not scale-free, 
thus ruling out a common theoretical simplification (e.g.. IChevalierlll982l ). 



3.1.2. Distance 



Comparing the angular size observed by VLBI with the radius inferred from optical, X- 
ray, and infrared observations gives a direct estimate of the distance to the remnant. Three 
estimates of the evolving shock veocity are currently available. 

First , the earliest optical spectra give initial velocities of 3500 ± 150kms~^ (based on 

Buill I2OO6I . and early spectra reported on VSNet), which we take as constant between the 
explosion (assumed to occur on Feb. 12.6 ± 0.15, or At = t — to = —0.23 ± 0.15 days with 
respect to tg, the time of optical maximum [Feb. 12.83]) and At = 1.47 days. 



-7- 



Second, shock velocities may be derived from the observed X-ray temperatures, using 
the standard strong shock approximation and assuming the standard mean molecular weight 
(/i = 0.6) for a fully ionized plasma. The X-ray emission is clearly complex, with indications 
of multiple temperatures and densities present at any given time. Here we are interested in 
the highest temperature gas, corresponding to the material heated by the decelerating shock 
shown in the radio images, and adopt single-temperature fits to the highest energy X-rays 
seen in a given observation. The X-ray temperatures at early times, when t he temperatures 



are h ighest, are best determined by the Rossi X-ray Timing Explorer (RXTE: ISokoloski et al. 



20061 ) . which responds most efficiently to higher-energy photons; at late times, as the shock 
slows and the corresponding temperature decr eases, the relative ly soft response of Swift 
provides the mo re relevant spec t ral co verage (JBode et al.l l2006bl ). We therefore use the 
measurements of ISokoloski et al.l (120061 ) before day 14.5, and those of iBode et al.l (l2006bl ) 
after day 13.5. The X-ray velocities are well fit by two independent power laws, as detailed 
in Tableland shown by the dashed lines in Fig. HI 



The third and final set of velocity estimates comes from infrared spectroscopy. iDas et al. 



(l2006bl ) measured the (deconvolved) full-widths at half maxima (FWHMs) and half-widths 
at zero intensity (HWZIs) of the observed line profiles of two infrared lines, O I and Pa/?, from 
1.16 to 47.04 days after the explosion. While the measurements for the two lines are in rough 
agreement, and the FWHMs agree with the optical and X-ray velocities discussed above, the 
FWHMs for both lines - and the velocities inferred from the optical and X-ray data - are 
systemat ically low (by a fac tor ~ 1.7) compared with the FWZI. This discrepancy was also 
noted by lEvans et al.l (120071 ). Fig. 2 of that paper shows that the difference between FWHM 
and HWZI results from the shape of the line profile: a strong, narrow, roughly Gaussian peak, 
surrounded by a faint, broad plateau. The narrow component suggests an exp anding wind 



or outfiow, while the broad plateau is more characte ristic of a thin shell (e.g., iLynch et al. 



20061 : iBertout fc MagnanI EqStI : Icill fc 0'Brlenlll999h . The consistent ratio of the FWHMs 
and HWZIs, and the roughly constant ratio of the peak fiux density in the two components, 
shows they remain closely associated for at least the first 40 days after the explosion. 

The early optical reports may well have missed a faint, broad line component. The much 
lower velocities derived from the X-ray temperatures, and their agreement with the FWHM 
of the narrow component of the infrared lines, is more puzzling. iTatischeff fc Hernanzl (120071 ) 
suggest that cosmic ray acceleration at the blast wave may explain both the velocity dis- 
crepancy and the relatively rapid deceleration of the shock. 

Here we assume the broad component of the infrared lines represents the shock velocity, 
as both the line shape and the extreme velocity suggest a thin, rapidly expanding shell. The 
scatter in the infrared measurements from day to day suggests some short-lived physical 



effects, or some uncertainty in the determinations of the hne widthsjj We therefore multiply 
the smooth fit to the X-ray and optical velocities (see above) by 1.7, to give a smooth curve 
which roughly reproduces the decline in the infrared HWZI. With these assumptions, the 
errors on the shock radii at the times of our VLBA observations are dominated by the error 
on the time of the explosion. 

Comparing these radii to the angular sizes corresponding to the outer part of the shock 
(6 = 17.5 ±0.5 mas at At = 20.8 days; 9 = 19.7 ±0.5 mas for At = 26.8 days) gives two fairly 
independent measurements of the distance: 



D20.8 
-D26.8 



'"20.8/^20.8 

'"26.8/^26.8 



2.392 ± 0.078 kpc 
2.492 ± 0.071 kpc 



where the error bars are statistical, and dominated by the uncertainties in 9. Taking the 
weighted average of these two distances gives a distance of 2.45 ± 0.05 kpc. 

The major uncertainty in this estimate lies in the relation between the sizes measured 
with VLBI, and the velocities inferred from other wavelengths. First, if the ring is a two- 
dimensional structure inclined at an angle i to the line-of-sight, that will increase the quoted 
distances by 1/ cosz: a factor 1.6 for i = 50°. Given the observed circularity of the ring, this 
is likely to be a < 10% effect. Second, we have assumed that the infrared velocities reflect 
those of the leading edge of the radio emission. Adding an additional 10% (rms) systematic 
uncertainty based on the possibility that the outer edge of the observed radio shock lies 
within the high-velocity region traced by the extreme infrared velocities, the final distance 
estimate is then 2.45 ± 0.05 ± 0.37 kpc, where the first error bars refers to the random and 
the second to the systematic errors. 



Previous distance estimates are summarized by iBodd (119871 ): see also Barry et al. (in 
prep.). These were based on spectral luminosity determinations, the presence or absence of 
absorption features at known distances, and a very rough rela ti on between t he distance and 
(cold) H I column density (see primarily iHjellming et al.lll986l ). iBodd (119871 ) concluded that 
the distance was likely to be around 1.6 ± 0.3 kpc. Our more direct estimate is higher but 
not enormously so. Further, since a significant inclination for the ring can only increase our 
distance estimate, the comparison implies again that the ring is nearly face-on. 



•^Fig. 2 of lDas et al.l (|2006bl ) shows error bars which appear sufficient to explain much of the scatter, but 
those error estimates are not given exphcitly in the text or in Table 2 of that paper. 
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3.1.3. Emission and Absorption 



Given the close agreement between the blast wave velocity evolution inferred from X-ray, 
optical, and infrared observations and the expansion of the radio ring, it seems worth checking 
whether a significant fraction of the radio flux density could b e thermal bremsstrahlung from 
the X-ray-emitting gas. Fits to the data of lBode et al.l (l2006bl ) (Table[2]) give the temperature 
of the post-shock gas as ~ 2.4 x 10'' K on day 20.8 and ~ 1.8 x 10'' K on day 26.8. These 
plasma temperatures are close to the observed peak surface brightness, which corresponds to 
an optically-thick black body of a few times 10^ K. However, a mass of more than 10~^ Mq 
would be required to make the thermal gas optically thick in the radio, and the corresponding 
thermal energy woul d be more than 10'^"^ ergs - an order of rnagnit ude above the total energy 



of the outburst (e.g., iHachisu fc Katdl200ll : ISokoloski et al.ll2006l ) - with a correspondingly 



high X-ray luminosity. Clumping of higher-temperature gas would reduce the required mass 
somewhat, but not enough to remove these difficulties. Thermal emission is thus unlikely to 
play a major role in the few-GHz radio flux at this stage of the outburst. 

The obvious alternative is synchrotron emission, which can easily explain the observed 
high surface brightness. Optically-thin synchrotron spectra generally range from z/~°'^ to 



,-1 



V ", providing a natural explanation for the ~ z/^"' spectrum of the western side of the 
ring. The eastern part of the ring, with a = —0.20 ± 0.02 {S^y oc z/'^) on day 20.8 and 
a = —0.27 ± 0.01 on day 26.8, is more puzzling. Foreground absorption in a 1/r^ medium 
would be expected to change rapidly over this time period, while the observed spectrum is 
fairly constant. Also, whereas any foreground absorption (or synchrotron self-absorption) 
should have caused the flux density to decrease with decr easing frequency below 1.7 GH z, in 
fact the 0.61 GHz flux density of 48 ± 2 mJy on day 20.3 ( Anupama fc Kanthariall2006l ) was 
higher than the 1.7 GHz flux density of 36.7 ± 0.3 mJy that we measured from the eastern 
part of the ring on day 20.8. Either the eastern side of the ring contains a complicated mix 
of emission and absorption, or the spectrum reflects a very unusual energy distribution of 
relativistic particles. 



Assuming the emission is synchrotro n - a conclusion shared by lO'Brien et al.l (J2006d ) 
- standard equipartition arguments (e.g.. lBurbidgdll959l ) suggest minimum energies in rel- 
ativistic particles and flelds of ~ 10"^^ (1 + /c) ergs, where k is the ratio of the energy in 
protons to t hat in electrons; /c ~ 1 00 is observed locally, while fc ~ 40 is predicted for strong 
shocks (e.g., iBeck &: Krausell2005l ). The equipartition magnetic fleld at the location of the 
peak flux density in the ring is o f order 0.03 (1 -|- k) G, comparable t o previous results 



based on much lower quality data (JBode fc Kahrull985l : iTaylor et al.lll989l ). The correspond- 
ing magnetic pressure is 2.7 x 10~^dyncm~^. As expected, these values are well below the 
total energy and pressure in the shock front. 
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3.2. The Eastern Component 

3.2.1. Emission Mechanism and Speed 

The steep spectral index (a ~ —0.67) and high surface brightness {Tb ~ 10^ K at 
1.7 GHz) of the eastern component indicate that the emission must be primarily synchrotron. 
The minimum (equipartition) energy is then ~ 0.15 x 10^^ (1 + fc) ergs, a factor 10 below 
that in the ring, while the corresponding equipartition magnetic field is comparable, ~ 
0.03 (1 + fc)^/^ G. 

If this component was ejected from the source at the core of the expanding ring, it 
is moving quite rapidly. Presumably it was ejected during or after the initial outburst, 
and is thus at most 27 days old when first detected by us (2006 Mar 11). At a distance 
of ~ 52 mas from the center of the ring, its projected velocity is > 2masday^^, which 
corresponds to > 8500/ sin z km s~^ (where i is the inclination to the line-of-sight, and taking 
a distance of 2.45 kpc). For an inclination i ~ 50°, the velocity of the eastern component 
is ~ ll,000kms~^. For comparison, the escape s peed from a white dwarf is ~ 6530(1 — 



M/1.456M0)-^/^kms-i flHansen fc Kawalei 



1994[). or - 12,600kms-i for M = 1.35 Mq 



Based on the contour plot of lO'Brien et al.l (l2006d ). the EVN detection on day 21.5 gives 



roughly the same speed. The eastern component had thus not decelerated significantly before 

day 27. 



3.2.2. First Appearance 



O'Brien et al.l (l2006d ) present a preliminary view of VLBA, EVN, and MERLIN ob- 
servations of RSOph in 2006, with the VLBI data covering days 13.8-28.7. Most aspects 
of our images agree reasonably well with their results. However, they detect the eastern 
component quite strongly (peak ~ 0.8mJybeam~^) at 5 GHz with the EVN on day 21.5, 
while we do not see that component with the VLBA on day 20.8. This apparent discrepancy 
is likely due to the more centrally concentrated uv-coverage of the EVN observations, which 
are correspondingly more sensitive to faint, extended emission. The eastern component was 
therefore already extended by day 20.8, as we directly observe on day 26.8. Further, our 
detection of the component at 5 GHz with the same instrument and observing time on day 
26.8 as our non-detection on day 20.8 shows that the eastern component was still brightening 
during this time. 

Our 1.7 GHz VLBA observations on day 20.8 do not show the eastern component. This 
result agrees with the non-detection of O'Brien et al. (2006) at the same frequency on day 
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20.5. Since the VLBA is much more sensitive to extended structures at this lower frequency 
than at 5 GHz (the angular response scales as A/i?, with A the wavelength, and B the length 
of the baseline under consideration), our non-detection at 1.7 GHz places strong constraints 
on the spectrum around day 21. Had the spectrum on day 21 been as steep as the spectrum 
we observed on day 27, we would have measured a 1.7 GHz flux density of > 1.7mJy, or 
> lOcr in our VLBA observations. Considering instead the 5a upper limit of 0.6 mJy /beam 
at 1.7 GHz, and comparing this limit to the EVN detection, we find that on day 21, the 
radio spectrum was actually either fiat or inverted (rising towards higher frequencies). The 
eastern-component spectrum therefore evolved significantly between day 21 and day 27. 



3.2.3. Absorption 

The spectral and flux density evolution, combined with the inferred large extent of the 
eastern component, suggest that external (free-free) absorption is important, as expected 
for a source embedded in a dense, ionized red giant wind. The primary challenge lies in 
explain ing the initial obscurat ion of the eastern component at a time when the ring is clearly 



visible. lO'Brien et al.l (l2006d ) invoke a special geometry to explain the visibility of the ring, 
as well as the late turn-on of the eastern component and the lack of any corresponding western 
feature during the first month. In their picture, the obscuring medium is the circumstellar 
material (assumed spherically symmetric) shed by the red giant companion. They assume 
the radio emission comes from a bright, slowly expanding equatorial ring and a bipolar 
shock structure. They take the approaching (eastern) side of the bipolar jet to be inclined 
about 30° to the line-of-sight, moving outward, and sweeping up absorbing material. In this 
picture, the eastern component of the jet becomes visible when the shock overtakes most of 
the absorbing circumstellar material. The receding jet remains behind more of the red giant 
wind, so that (western) component is not expected to be seen until later. They suggest that 
the more distant (eastern) part of the ring is visible almost instantly because the jet has 
swept away the densest, high-obscuration part of the red giant wind along that line-of-sight. 

One can make this more quantitative by requiring that the free-free opacity r due to the 
red giant wind be > 1 at the times and frequencies where the eastern component is not seen, 
and ^ 1 at the times and frequencies where it is. The opacity predicted for a fully ionized, 
spherically symmetric, 1/r^ wind are shown as a function of projected distance from the 
center in Fig.O together with arrows showing the rough constraints from the 1.7 and 5 GHz 
VLBI observations. This simple model works remarkably well. Assuming one proton per 
electron and an electron temperature of 10^ K, and taking the wind velocity v^ = 20kms~^, 
the inferred mass loss rate is ~ 1 x 10^^ MqJt^^ for an inclination i ~ 50° for the eastern 
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component's velocity vec tor. This mass loss rate is qu ite reasonable, and within a factor 
of two of that derived by lO'Brien. Bode. &: KahnI ( 19921 ) (after correcting for the difference 
in the assumed distance) based on shock models fit to the X-ray emission from the 1985 
outburst. 

Absorption can thus easily explain the late appearan c e of th e eastern component; but 
the central ring is more challenging. In the lO'Brien et al.l (J2006d ) model, the observed ring 
should represent the projection of a highly- inclined, two-dimensional circle, i.e., an ellipse. 
Our data show that the ring in fact is quite circular (Fig.[2]), making explanations involving 
strong deviations from spherical symmetry difficult to m aintain. The flat spec trum of the 
eastern part of the ring is also not a natural result of the lO'Brien et al.l (l2006d ) model. 



Another possible explanation for the abrupt appearance of the eastern component is 
that it was intrinsically faint until around day 21. The eastern component could represent a 
terminal shock of the kind associated with the ejecta of many active galactic nuclei, micro- 
quasars, and protostars (in the form of Herbig Haro objects), although this explanation has 
some difficulties since for RS Oph, the circum-binary density falls off rapidly with distance 
rather than abruptly increasing. 

The ejecta could also have brightened upon moving into the relative vacuum beyond 
the outer edge of the wind nebula. As in the current outburst, the expanding shock in 1985 
swept up the CSM, leaving behind a relative vacuum. Over the past 21 y ears, the red giant 
has been re-filling this vacuum with a wind moving out at ~ 20kms~^ (jWallersteiru Il958 



Dufay et al.l Il964l : iGorbatskiiil Il972l ). The collimated ejecta would have reached the edge 
of this re-filled region in ~ 14days (fu;/20kms~'^) (fjet/H, OOOkms"^) (where Vjet is the 
velocity of the collimated outfiow), around the same time that it brightened in the radio. If 
the red giant wind can carry outward and compress (and hence strengthen) the magnetic 
field at large radii, relativistic particles might produce copious synchrotron emission upon 
encountering the compressed field. Alternatively, if the relativistic plasma was initially highly 
confined by the circumstellar gas (and thus synchrotron self-absorbed), it might expand and 
become visible at lower radio frequencies upon reaching the edge of the confining mediumo 
One might then expect the eastern component to become relatively stationary at late times 
rather than continuing to move out rapidly. 



*Note that the azimuthal extent of the eastern component requnes a fairly wide-angle 'jet'. If the jet 
were intrinsically narrow until reaching the edge of the red giant wind, the azimuthal velocity would have to 
be > 17, 000 km s""'^ to match the observed azimuthal extent on day 21. 
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4. Conclusions 

Dual-frequency VLBA observations made 20.8 and 26.8 days after the February 2006 
outburst of the recurrent nova RS Oph have shown that the bulk of the radio emission at 
this stage comes from an expanding synchrotron shell. The shell's increasing angular size 
matches the shock velocities inferred from X-ray observations and measured through infrared 
spectroscopy, confirming the basic picture of a decelerating shock between the material flung 
off in the explosion from the surface of the white dwarf and the circumstellar material 
produced by the wind of its red giant companion. This gives a direct distance estimate of 
~ 2.45 ± 0.05 ± 0.37 kpc, the most accurate yet determined for this source. 

This evolving shock is far from simple. The shell evolution is not self-similar, with the 
shell thickness increasing from ~ 15% to ~ 26% between the two epochs. In addition, the 
brightness distribution around the shell is highly asymmetric, with the eastern side many 
times brighter than the western side. The eastern side of the ring also has a much flatter 
radio spectrum. To our knowledge, such a spectrum has been seen only once before for a 
high-brightness resolved source, in the radio images of the e xpanding shell associated wit h 



the peculiar X-ray outburst of the emission-line star CI Cam (JMioduszewski fc Ruperul2004l ). 
Finally, the flux density of the western part of the shell is increasing with time, in contrast 
to the optically-thi n behavior of the ej ecta-driven shocks in physically similar systems such 



as supernovae (e.g. 



Weiler et al.ll2002[ ). 



Our VLBA images also reveal an additional eastern component that is well separated 
from the dominant ring. 'Jets' (generally just elongations or separate components in barely- 
resolved sources) have been seen in previous outbursts of this and other symbiotic stars. The 
VLBI observations of RS Oph constitute the flrst unambiguous evidence for a synchrotron 
jet emanating from a white dwarf system. The inferred speed is close to the escape speed 
from a high-mass white dwarf, in accord with the generally good agreement between directed 
outflow and escape speeds for physical systems ranging from low-mass protostars to black 
hole binaries and active galactic nucleo Similarly, synchrotron emission is characteristic of 
fast jets, where emission far from the powering source is generally assumed to derive from 
an interaction between the flow and the surrounding medium. The association of a highly 
directed and a quasi-spherical outflow is more unusual, but certainly not unprecedented; 
for instance, the jet model for gamma-ray bu rsts, and their association with core-collapse 



MacFadyen fc Woosleylll999l ). has brought the physical connection between 



supernovae (e.g. 

quasi-spherical and highly directed outbursts into much greater prominence. 



^In supernovae and gamma-ray bursts, the relevant escape speeds are those of the compact objects which 
are formed, not those of the progenitor stars. 
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The eastern component that was so prominent on day 27 appeared to "turn on" around 
day 21. Although this turn-on could be due to decreasing absorption by the circumstellar 
medium as the eastern component moved outwards, there are several difficulties with this 
picture. Intrinsic brightening, possibly associated with the component reaching the edge of 
the region re-ffiled by the red giant wind since the last nova outburst in 1985, may also play 
a role in the turn-on. 

The preferred east-west axis for the observed asymmetries is consistent with VLBI ob- 
servations on day 7 7 of the last (1985) outburst, which showed a strong east /west extension 



(iTaylor et al.lll989l ). The persistent east- west asymmetry suggests that this preferred axis 
reflects an intrinsic asymmetry in the binary system, rather than (for instance) some sta- 
tistical irregularity in a particular explosive event. The binary orbit provides one possible 
source of the asymmetry. Although a highly anisotropic density structure around the binary 
cou ld coUimate the ou tflow, with the 'jet' naturally focused along the axis of least resistance 



(cf. iLloyd et al.lll993l ). one would then expect the observed ring to correspond to the high- 
density 'waist' of this structure. The ring is in fact circular in projection and brightest along 
the trajectory of the eastern component. Whatever the coUimation mechanism in RS Oph, 
we may for the ffist time be seeing the formation of a white-dwarf jet, as well as the most 
compelling demonstration to date of both its high initial speed and the physical emission 
mechanism (synchrotron emission). 

We are grateful to the VLBA scheduling committee and the array operations staff for 
the efforts which made these observations possible. As with most research on bright variable 
stars, this work has greatly benefited from the regular monitoring (and quick reporting) 
of both amateur and professional variable star observers, generally publicized through the 
American Association of Variable Star Observers (AAVSO) and the Variable Star Network 
(VSNet). We are further grateful to the referee, Mike Bode, for a number of useful comments, 
as well as for organizing an exceptionally productive workshop on this source in the summer of 
2007. J.L.S. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship 
under award AST 03-02055. This research has made use of NASA's Astrophysics Data 
System; the SIMBAD database, operated at CDS, Strasbourg, France; and the USNOFS 
Image and Catalogue Archive, operated by the United States Naval Observatory, Flagstaff 



Station (http: //www. nofs .navy .mil/data/fchpix/] 



Facilities: VLBA (). 
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Fig. 1. — VLBA total intensity images of RS Oph, from 2006 March 5 (day 20.8; upper) 
and 2006 March 11 (day 26.8; lower), (a) 5.0GHz images (left column), convolved to a 
resolution of 7.6 x 7.6 mas. Contours are ±0.25 mJybeam"^ x 2"/^, n = 2, 3, . . ., and the color 
scale ranges from to 6mJybeam~^. The peak flux densities in the two images are 5.9 ± 
0.25mJybeam^^ on day 20.8, and 4.3±0.11mJybeam~^ on day 26.8. Note that the 5.0 GHz 
images are inherently less sensitive to large structures, due to the difference in interferometric 
baseline coverage as measured in wavelengths. The total 5.0 GHz flux densities, uncorrected 
for this effect, are 41.3 and 26.3 mJy. (b) 1.7 GHz images (right column), convolved to a 



resolution of 13 x 13 mas. Contours are ±0.12mJybeam 



-1 



2n/2 



n 



2,3, 



and the 



color scale ranges from to 13mJybeam~^. The peak flux densities in the two images are 
13.0 ± 0.12 mJybeam"^ on day 20.8, and 10.3 ± 0.13 mJybeam"^ on day 26.8, with the error 
bar reflecting the off-source rms noise. The total flux densities are 43.1 and 30.1 mJy. The 
cross to the northeast in each panel marks the optical position of RS Oph, from NOMAD 
fsee 512.11): the size of the cross represents the NOMAD error bars. 
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Fig. 2. — Tilted circular rings superposed on the 5.0 GHz image from day 20.8, for inclinations 
of 0° (face-on) and 50°, using radii of 31 and 29 milliarcseconds, respectiv ely. The orbital 



inclination is thought to be ~ 50° fJQuiroga et al.ll2003l : iBrandi et al 



20071 ). The color scale 



for the image ranges from 0.25 (black) to 3.865 (red) mJybeam ^. 
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Fig. 3. — The radius of the ring as a function of time since the optical peak on 12.83 Feb 
2006 UT. The soh d circles represent the outer radii of the fit annuli (see text). The triangles 
shows the data of lO'Brien et al.l (l2006d ). assuming a 1 mas error bar (see text); these size 
estimates are only approximate, and the disagreement between the t wo groups is lik e ly due to 
differences in the methods of their derivations (e.g., it is likely that lO'Brien et al.l ( l2006d )'s 
values refer to the mid-point rather than the outer edge of the ring). The dashed line shows 
the expansion detailed in Table[2|, scaled up by a factor 1.7 to match the width of the broad 
infrared spectral component (see § l3.1.2p . 
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Fig. 4. — The velocity of the shock vs. time, as derived from infrared, optical, and X-ray 
data. The observed widths of the O I line are shown in green, and those of the Pa/3 line in 
light blue (cyan); the open squares represent the full- widths at half-maximum, and the solid 
squares show the half-widths at zero intensity (H WZIs). The blue dots and error bars are 
shock velocities inferred from RXTE'K.-isy da ta inlSokoloski et al.l (120061 ). while the red ones 
are similarly derived from Swift X-ray data in lBode et al.l (l2006bl ). The blue and red dashed 
curve s represent pow e r-law fits (see Table |2]) to the observed X-r ay velocities , ignorin g data 
from ISokoloski et al.l (120061 ) taken after day 14.5, and data from iBode et al.l (J2006bl ) taken 
before day 13.5. The dashed black horizontal line shows the velocity for the earliest part of 
the outburst, taken from early optical reports. As discussed in § 13.1.21 the solid black, blue, 
and red curves show the result of multiplying the X-ray and optical fits by 1.7, to match the 
infrared HWZIs. These solid curves represent our best estimate of the shock velocity. The 
vertical dashed lines show the dates of our VLBA observations. 
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Fig. 5. — The free-free opacity predicted for a 1/r^ wind from the red giant, at 5 GHz (blue) 
and 1.7 GHz (red), as a function of the (observed) projected distance from the red giant. 
The red and blue arrows show the rough constraints from the VLBI imaging of the eastern 
component — the opacity must be < 1 when that component is not detected, and > 1 when 
it is detected. The horizontal placement of these limits assumes D = 2.45 kpc; a smaller 
distance would shift these to the left (implying a lower M/vw), while a larger distance would 
shift these to the right (implying a higher M/vw). The solid curves show a rough fit to these 
limits, corresponding to M/v^, = 1 x lO~^M0yr^^/2Okms~^ for an inclination of i ~ 50°. 
The dotted lines show the opacities for a 20% lower M/vw, while the dashed lines show the 
opacities for a 20% higher value. These calculations take the electron temperature to be 
10^ K, and assume one proton per electron. 
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Table 1. 



Epoch 



Feature 



Si. 7^ 

(mjy beam^"'^) 



S, 



corr, 5.0 



(mJy beam 



-1^ 



a 



2006 Mar 05 



2006 Mar 11 



Entire ring 
Ring-E<^ 
Ring-W^ 

Entire ring 

Ring-E'i 

Ring-W^ 

E.C.'i 



43.1 ±0.3 

36.7 ±0.3 

6.4 ±0.3 

41.3 ±0.2 

33.3 ±0.2 

8.0 ±0.2 

4.2 ±0.2 



32.2 ±0.6 
29.5 ±0.6 
2.8 ±0.6 
28.8 ±0.2 
24.8 ±0.2 

3.8 ±0.2 

1.9 ±0.2 



-0.27 ±0.02 
-0.20 ±0.02 
-0.76 ±0.20 
-0.33 ±0.01 
-0.27 ±0.01 
-0.67 ±0.01 
-0.71 ±0.10 



^Flux density at 1.7 GHz. 

^Flux density at 5.0 GHz, corrected for flux density missing as found by 
simulating 5.0 GHz observations of the 1.7 GHz model (see text). The observed 
flux densities have been muhiphed by 1.07 on 2006 March 05, 1.10 for the ring 
on 2006 March 11, and 1.21 for the eastern component on 2006 March 11. 

'^Spectral index a, with the flux density going as the frequency to this power 
{S, oc z/"). 

■^Ring-E refers to the eastern (bright) part, and Ring-W to the western (dim) 
part, of the ring. E.G. refers to the eastern component. 
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Table 2. 



Time coverec 


i by fit^ 

t2 

(days) 




Fit parameters'^ 




(days) 


(A.U.) 


«d 


xVN,,,j,^ 


-0.23 ±0.15 


1.47 


3.44 ±0.34 


= 0.0 




1.47 
14.0 
14.0 


14.0 
20.8 
26.8 


15.49 ±0.17 

cr 7r|+0.039 
^■"-'-0.038 

q qc+0.063 

y-yD-0.056 


-0.350 ±0.018 

-n c;94+0.022 

U.0Z4_QQj^g 

n c;9A+0.022 


2.94/3 

2.58/3 
2.58/3 


J- f 

'^expl 

J- f 
'^expl 


20.8 
26.8 


24.63 ±0.38 

28.89+[J:| 







'^Times ti to ^2 over which the fit holds, measured in days with respect 
to the time of optical maximum to =Feb. 12.83. 

'^Results of least-squares power-law fits to the X-ray data, parameterized 
in terms of the distance traversed by the shock i?2i, and the power-law 
index a. Note that the distances have not yet been scaled to match the 
broad infrared spectral component - see § 13.1.21 

■^The distance the shock has traveled between ti and t2, in astronomical 
units (A.U.). Note that these distances have not yet been scaled to match the 
broad infrared spectral component - see § 13.1.21 Error bars are la based on 
standard propagation of errors, and are primarily statistical. The exception 
is the uncertainty in the explosion date {texpi = —0.23 ± 0.15), which is a 
best guess as to the likely range of dates. 

'^The power-law index a of a fit to the velocities: f oc t". Error bars are 
la. 

^X^ and the number of degrees-of-freedom for the fit. Note that x^ has 
not been normalized — X^ ~ Nd.o.f. indicates a good fit. 

^These rows give the total distance traveled between the (uncertain) date 
of the explosion texp^ and the dates of our VLBA observations (20.8 and 
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26.8 days). The error bars are la based on a simple propagation of the error 
bars for each segment of the fit, assuming the error bars on the individual 
segments are independent. Note that these distances have not yet been 
scaled to match the broad infrared spectral component - see § I3.1.2[ 

Note. — The results given here have not yet been scaled to match the 
velocities measured for the broad infrared spectral component - see § I3.1.2[ 
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